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We examine the dynamics of bipartite entanglement between a two-level atom and the electro- 
magnetic field. We treat the Jaynes-Cummings model with a single field mode and examine in 
detail the exact time evolution of entanglement, including cases where the atomic state is initially 
mixed and the atomic transition is detuned from resonance. We then explore the effects of other 
nearby modes by calculating the exact time evolution of entanglement in more complex systems 
with two, three, and five field modes. For these cases we can obtain exact solutions which include 
, the strong-coupling regimes. Finally, we consider the entanglement of a two-level atom with the 

' infinite collection of modes present in the intracavity field of a Fabre-Perot cavity. In contrast to 

^SJ ' the usual treatment of atom-field interactions with a continuum of modes using the Born-Markov 

approximation, our treatment in all cases describes the full non-Markovian dynamics of the atomic 
O ' subsystem. Only when an analytic expression for the infinite mode case is desired do we need to 

' make a weak coupling assumption which at long times approximates Markovian dynamics. 
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INTRODUCTION 



. In recent years the quest for a quantum computer has generated much interest in quantum information processing 
' [H and the need for a deeper understanding of the distinguishing quantum features of many body systems, such as 
^ ^. quantum coherence and entanglement. Recent progress in many areas of physics, especially atomic-optical, condensed 
matter, low temperature and mesoscopic physics, has been brought about by new techniques (in manipulating atoms 
and light) and high precision measurements and fabrication of new materials. The drive for better preparation, manip- 
CN ■ ulation and control on the experimental side is matched by an adjustment of focus and emphasis of theoretical studies. 
K*" ' There is a trend to go back to simpler models (e.g., Bose-Hubbard for BEC) but explore the parameter regimes and 
. physics hitherto unknown or unclear (such as the BCS-BEC transition, BEC-Mott transition) but made possible only 
' recently by the high precision and highly manipulative experiments. With these efforts new demands on experimental 
designs (such as accurately positioning an atom in a cavity 0, 01) for new devices (such as single atom laser 0, 
and new theoretical issues (such as quantum fluctuations and correlations, decoherence and disentanglement) come 
to the fore. 

It is in this spirit that we explore some of these newly introduced issues in atom-field interaction, with the stated 
] purpose of better understanding the effects of non-Markovian dynamics and strong coupHng between a two-level 
atom (2LA) and the electromagnetic field (EMF). In this paper we will seek first to give a detailed, general, and exact 
treatment of the dynamics of entanglement in the Jaynes-Cummings model. From there, we expand our considerations 
to include the effects of other, nearby field modes, which can, in principle, add new features to the evolution of 
^ ■ entanglement when coupling is strong enough. Finally, we will explore interaction between the atom and the full, 
5^ ', infinite collection of modes in the intracavity field. In all these cases we solve for the exact dynamics j2l| keeping the 
full unitary dynamics of the composite system, so that our treatment can capture the non-Markovian behavior of the 
atom (considered as a subsystem) even in the infinite-mode case. In this approach, then, we wish to give a detailed 
description of the much-studied, single-field-mode Jaynes-Cummings model and connect this to strong-coupHng and 
non-Markovian phenomena in more general models. We want to adopt an approach which can best preserve the 
quantum coherence and entanglement of the system and include the full interplay of the subsystems involved (or 
back-action from its environment, if any one such subsystem merits special attention), treated self-consistently. The 
infiuence functional approach we take here is particularly suitable for such requirements, and they are necessary 
for a correct treatment of disentanglement, because in general it is an even more intricate and delicate issue than 
decoherence. 
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Simpler atom-field models, such as that of a 2LA interacting with one single EMF mode described by the Jaynes- 
Cummings model, have solutions in closed form. Such closed form solutions offer a good comparison with results 
involving approximations. Quantum entanglement shows up as 'collapse and revival' phenomena, and has been 
studied in depth by Pheonix and Knight 0, Q and Gea-Banacloche 15 years ago for initial pure states and Bose et 
al [§| for initial mixed states. We analyze this problem to obtain the exact time evolution of the entanglement (rather 
than just a bound on or average of entanglement) for general, mixed initial atomic states, and we explore the effect of 
detuning between the atomic transition frequency and the cavity mode frequency. This gives us a detailed and exact 
picture of the dynamics of entanglement in the well-known Jaynes-Cummings model. 

Moving on to more complex models, we study the case of a 2LA interacting with two, three, or five EMF modes. 
We present an exact solution of the time evolution of entanglement showing the effect of these additional models for 
an arbitrary, pure initial atomic state. Using a perturbative approach we give leading order corrections to Jaynes- 
Cummings dynamics from other, nearby field modes in the case of strong atom-field coupling. There is a recent 
proposal to use the 2LA as a control qubit and the two EMF modes as a target qubit in the making of a CNOT gate 
Our calculation of entanglement between an atom with an initial mixed state and multiple field modes may be 
useful for these proposed studies. 

In support of the new needs mentioned above, our work aims at probing the strong-coupling regime. Strong coupling 
is defined as g >> k, 7, where k is the cavity damping rate and 7 is the rate of spontaneous emission from the atom. 
In this regime the coherent evolution of the atomic and intracavity field state is fast compared to the dissipation rates 
in the system, so coherent quantum effects become important This is the regime we explore in treating the few 
mode cases. 

Beyond the strong-coupling regime, there is also the regime of super-strong coupHng g ~ A, where A denotes 
free spectral range of the resonator in angular frequency. While, conventionally, strong coupling is achieved through 
recycling of the light by means of a high Q cavity, in the super strong regime the coupling between atoms and light is 
already strong during one round trip in the resonator. Thus the atom can affect the mode structure, in addition to the 
occupation of the modes, as is the case in conventional Cavity QED. (This, in our language is when the back-action 
of the atom on the field begins to become important.) According to this super-strong-coupHng regime can be 
achieved for a modest number of ultracold atoms in an optical lattice formed by standing waves in an optical resonator. 
Our studies of cases with a few field modes cover this super-strong regime for the case of a single atom, showing new 
behavior in this regime, and they apply to certain situations for the multi-atom scenario which is experimentally 
realistic. 

Finally, we study the interaction of a 2LA with a cavity and its full, infinite collection of EMF modes in certain 
limits, going beyond a simple one- or few-mode treatment. We do not need to invoke the Born-Markov approximation 
and our approach generally appHes to strong-coupHng regimes. Our calculations yield formal solutions for the evolution 
of entanglement in terms of integrals, which admit numerical solutions in the general case, but in order to obtain 
simple analytic solutions we explore a solution for late-time behavior that is perturbative in the atom-field coupling, 
and, so, those expressions are limited to the weak-coupling regime of cavity QED and exhibit behavior corresponding 
to Markovian expectations. 

The second aim for our work is to probe into the non-Markovian regimes. Since our system, comprised of one atom 
together with the modes of the EMF, is closed, the dynamics of the complete system is unitary. One may wonder 
where non-Markovian dynamics arises. This is a reasonable question because non-Markovian behavior is usually 
associated with open system dynamics, which is the case when one focuses on one of the subsystems, often the 2LA, 
and introduces some coarse-graining in the other subsystem, such as EMF, whence one begins to see dissipative and 
stochastic behavior in the open system. However, strictly speaking, non-Markovian just refers to history-dependent 
processes (involving memories), and it is not just for open systems, in the following sense: It is common knowledge in 
non-equilibrium statistical mechanics that for two interacting subsystems the two ordinary differential equations 
governing each subsystem can be written as an integro-differential equation governing one such subsystem, thus 
rendering its dynamics non-Markovian, with the memory of the other subsystem's dynamics registered in the non- 
local kernels. All these are happening in a closed system except now one has shifted the attention to one of its 
subsystems, [l^l When an exact solution for the combined system is available, as in many of the cases studied here, 
one can derive the non-Markovian dynamics in the corresponding cases when the above mentioned procedures are 
taken, or if the non-Markovian dynamics is available from an open system consideration, one can identify the coarse 
graining measures introduced. The commonly invoked Markovian approximation gives rise to dynamics farther remote 
from (meaning, less accurate than) the exact solutions. 

We must add that investigations into the strong-coupling and non-Markovian regimes are two rather difficult tasks, 
even for cavity QED conditions. This study takes only a small step toward them. Our aim is to gain some intuitive 
understanding of the evolution of quantum entanglement in these new territories with a familiar and relatively solvable 
model. Much more work is needed to understand the general features for more complicated systems in these regimes. 

In Sect. II we specify the details of the physical system we wish to treat and discuss how we will quantify the 
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entanglement present. Sect. Ill discusses the method we adapted to find the evolution of the system and the effects 
of initial conditions assumed. We then treat the specific cases of different collections of EMF modes in Sect. IV, 
including a field consisting of a single mode, a few modes, and the full collection of modes present in an ideaHzed 
optical cavity. We end in Sec. V with a discussion and a summary. 

II. SYSTEM AND QUANTITIES UNDER INVESTIGATION 

A. Hamiltonian and Initial States Considered 

The system we consider is a two-level atom (2LA) coupled to M electromagnetic field (EMF) modes. We adopt the 
dipole and rotating-wave approximations, neglecting motion of the atomic center of mass. The Hamiltonian of this 
system, as given by is 

M 

H = hiooS+S- + ^ TiWkPh + h{gkS+b + glS-b'^) = Hq + Hi (1) 
fc=i 

where 

M 

Ho = huoS+S- + hjob'<b (2) 

is the free Hamiltonian and 

M 

Hi = Y. ^^kP'^ + K9kS+b + glS-b^) (3) 
fe=i 

is the interaction Hamiltonian. Here is the (bare) frequency of the atomic transition, lj^ is the frequency of the /c*'' 

field mode, gt is the coupHng of the fc*'' mode to the atom, and 5k = tok — ^o- Since gk oc l/y/uJk implies gj — gk^J^, 

we may choose to define things such that all values gk are real. Where the atom and the electromagnetic field modes 
in question form a closed system, the dynamics of the total system is unitary. For the M = 1 case, this is just the 
famihar Jaynes-Cummings model with no dissipation. 

We will assume a separable initial state described by the density matrix 

x(0)=/5a(0)®|0)(0| (4) 

where all electromagnetic field modes are in the vacuum state |0) and the atom may be in an arbitrary (possibly 
mixed) state, described by the atomic reduced density matrix pa{0)- If Pa{0) is pure then x{t) is pure for all times 
due to the unitary evolution. 

Tracing over the field, the reduced density matrix of the 2LA at time t is given by pA{t) = TrF{x{t))- As is well 
known, it can be represented by a point in the Bloch sphere. We will use the spherical polar parametrization of the 
sphere, so that pA can be specified by the triple {r,9,(l)), and, as usual, pure states lie on the surface of the sphere 
with r = 1. We choose (1,0,0) to represent the excited state of the atom |e) and (l,7r, (/)) to represent the ground 
state \g). 

In this work we seek to quantify the bipartite entanglement between the atom and the electromagnetic field (con- 
sidered as a whole). We will not look at entanglement with individual modes of the field separately, nor will we 
explore entanglement among the modes of the field. A local unitary operation of the form U = tlAtom ^ U Field does 
not change the amount of entanglement in the system and the time evolution due to Hq, which relates the 
Schrodinger picture to the interaction picture, is such a local unitary operation; therefore, we may work entirely in 
the interaction picture and compute the entanglement of the interaction picture state directly, as though it were the 
Schrodinger picture state. 

B. Measures of Entanglement 

When deahng with pure states, the Entropy of Entanglement (EOE) quantifies the amount of entanglement [l5| . 
The EOE is defined to be the von Neumann entropy of the reduced density matrix of the system. 

^vn{x) = S{TrF[x]) ^ S{pa) = - ^Pj ^og^ipj) (5) 
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where the pj represent the eigenvalues of the reduced density matrix pA- The EOE does not depend on which 
subsystem is traced over. We will always trace over the field degrees of freedom. 

When the total state of the system is mixed, the EOE can no longer be used as an entanglement measure. While 
many measures of entanglement have been defined, including the Entanglement of Formation, the Distillable En- 
tanglement, the Entanglement Cost, the I-Concurrence, and others, most of these involve an infimum or supremum 
over a large set (often the set of projector decompositions of x) that is in general difficult or impossible to evaluate 
analytically and quite costly to compute numerically pJi]. We choose to focus on an entanglement monotone known 
as the Logarithmic Ne gati vity (LN) jT^ because we can give explicit expressions for the value, and it can be easily 
calculated numerically |23l |. 

The LN is defined as 

£Ar(x) = log2||x^"||i (6) 
where is the partial transpose in the basis B of the density matrix and 



O =Tr{VOOU (7) 



is the trace norm of the operator O. The partial transpose in the basis B = p} is defined such that if 

X= I] X3,k,l,m\j){l\A<»\k){m\F (8) 
j^k,l,m 

then 

X^^ = X3,k,i,m\j){l\A<E)\m){k\F (9) 

j.k,l,rn 

SO that the second subsystem has been transposed. While the operator produced by this procedure will depend on 
which subsystem is transposed and in which basis, the trace norm depends only on the eigenvalues, which will be 
independent of these choices • 

In the case that the total state of the system is pure, the amount of entanglement can be determined from the 
reduced density matrix for one of the subsystems, as is clear from the definition of the EOE. In this case one may 
show, using the Schmidt decomposition, that the LN becomes 



£M^iog, VVpJ (10) 




again with pj as the eigenvalues of the reduced density matrix. When the total state of the system is mixed, however, 
knowing only the reduced density matrix is not sufficient (since, for example, the reduced density matrix of a maximally 
entangled state is the same as when the total system is in a separable, completely-mixed state). Thus, if we wish to 
track the evolution of entanglement for a total mixed state, we must keep more information than just the reduced 
density matrix. 



III. ENTANGLEMENT EVOLUTION 



A. Methods of Solution 



For M=l we have the Jaynes-Cummings model, and one can exactly diagonalize Hj in the basis known as the dressed 
states. The evolution of x in this basis is quite simple, and knowing x{t) we can then calculate the entanglement 
using the LN in general and the EOE if x(0) is a pure state. For more than one EMF mode, one can still diagonalize 
the Hamiltonian in the same way. Our initial state has all the electromagnetic modes in their vacuum, so the energy 
arises from the atomic state. This confines the initial state to a subspace spanned by the zero- and one-excitation 
energy eigenstates of Hq. Because the total system evolves via a time-independent Hamiltonian and [Hq, Hi] = 0, we 
know that the system will remain in this subspace. We then only need to diagonalize Hj in the subspace of degenerate 
one-excitation eigenstates of Hq. When there are M field modes, this subspace has dimension M + 1. While this 
diagonalization is straightforward numerically, doing this analytically requires the solution to a polynomial of degree 
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M + 1, so for M > 3 there is in general no closed form solution. In any case, having obtained the eigenvalues the 
solution is otherwise exact. 

For a pure initial state x(0), however, we only need the reduced density matrix pA to find the entanglement. 
Anastopoulos and Hu [l3| found the solution for the evolution of pA{t) in these types of systems with these initial 
conditions. Given an initial atomic state (1,6*, 0) at time t the reduced density matrix of the atom interacting with 
M modes is given by 



where 



PA{i) 



u{t) = 



^u*u{l + cos{0)) iue*"* sin(6l) 
iM*e-''^sin(6') 1 - + cos(6')) 



Cr^ represents the inverse Laplace transform, and 



1 



c+ioc 



-' Z + iOJi 

The entanglement between the atom and the field modes depends only on the eigenvalues of pA^fy, 




so it depends only on the norm of u(i). For M field modes 



\u(i)\ 



M ^ 



E 

J=0 



with the Zj values being all the solutions to the equation 



M 



M 



M 



l^k 



0. 



fc=l fc=l 

Using Eq. (fT4l) we can compute the FOE or the LN. The LN takes the particularly simple form 



£m = log2 ( 1 + 2cos^ ( - ) Vkp- 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 
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B. Dependence on Initial Conditions 



In examining the evolution of entanglement in this system, one of the primary questions is how it depends on initial 
conditions. Given the Hamiltonian and the class of initial states we are considering, the entanglement at any point 
in time is completely independent of the parameter 4> on the Bloch sphere for the initial state. While this is clear 
for pure states from the form of Eq. lfT4|) . we show this is true for mixed initial states as well. The unitary operator 
U^{a) = e-'-f^oa/(;ia;o) acting on an initial state of the form Eq. ^ will simply shift the angle 4> of the atomic state 
by an amount a. Consider two initial states, x(0) and x'(0) = U(f,{a)x{Q)Uij,{a)\ which differ only in their coordinate 
(j) on the Bloch sphere. Since commutes with Hi, we know 



(18) 



Because U,p{a) represents a local operation, we also know that the degree of entanglement in the states x{t) and x'{t) 
must be the same. 

For the case where the initial state of the system is a pure state, we can additionally conclude that, while the value 
of the entanglement at a given time depends on the 9 coordinate of the initial state, qualitative features like the times 
at which local maxima and minima of entanglement occur will be the same for all such initial pure states. The easiest 
way to see this is by looking at Eq. lfT7|). since the value of \u{t)\ at which the LN reaches extrema does not depend 

on e. 



(a) r = 1 



(b) r = 1/2 



(c) r = 



Figure 1: LN as a function of time t and the initial state value of on the Bloch sphere for states with different values of r, 
where g is the value of the atom-field coupling. In the last case, of course, all values of 8 correspond to the same state. 



EoE 




Figure 2: The dependence of entanglement on the parameters of the Hamiltonian, the detuning 5 and atom-field coupling g. 
The odd numbered minima in the entanglement (as a function of time t) increase and eventually disappear as the detuning 
increases. 

IV. SPECIFIC CASE STUDY 

A. Single Mode — The Jaynes-Cummings Model 

We begin with the case of a single mode of the field (M = 1) interacting with the atom in the familiar Jaynes- 
Cummings model. The Hamiltonian now depends only on the atom-field coupling g, the field mode frequency uj, 
and the atomic transition frequency lu + S. First, consider a system in which the atom is initially in a pure, excited 
state and the field mode is resonant with the atomic transition. In this case, the system undergoes Rabi oscillations 
between the separable states |e)|0) and passing through a maximally entangled state each way. As a result, 

the entanglement has the relatively simple oscillatory behavior with a time scale set by g. For an initial pure state 
we have 

u{t) = e-'("+i)* 

where Aa = \/ 5^ -\- Ag"^. For initially mixed atomic states, the fully density matrix for the bipartite system may be 
computed using the dressed states, as discussed in Sect. IIII Al The variation in the behavior of entanglement in the 
6 = case for different initial conditions is shown in Fig. [TJ As stated above it does not depend on 0, and many of 
the qualitative features such as the times at which maximal and minimal entanglement occur do not depend on 9. 

The other question is how the entanglement generated by the dynamics of the system depends on the parameters 
of the Hamiltonian. As is well known for the Jaynes-Cummings model, the time scale of the evolution is the vacuum 
effective Rabi frequency ^^/W+4g^ . What's more interesting is the dependence on the detuning 6 shown in Fig. [H 



Aa + 5 
2Aa 



Aa- 6 
2Aa 



(19) 
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Figure 3: The Bloch sphere for the one-excitation manifold Ho. The vector representing the quantum state of the system 
lip) precesses in a circle about the vector for the Hamiltonian Hi as time advances. When the state vector encounters either 
pole of the sphere, the state becomes separable. If the state vector crosses the equator of the sphere, the state is maximally 
entangled at that instant. The ratio of the detuning 5 to the atom-field coupling g determines the angle the Hamiltonian vector 
makes with the vertical axis. If 5 = the Hamiltonian vector is horizontal, and as 5/g becomes large the Hamiltonian vector 
approaches the vertical. The vector for the state starts on the pole of the sphere (due to the initial condition) and travels 
around a circle centered on the Hamiltonian vector; thus, the circle of precession is a great circle when 5 = 0, and the circle 
becomes very small when S/g is large. 



For an initial pure state, we can understand the time evolution of entanglement by focusing on the two level system 
composed of the one-excitation states of Hq. This two level system can again be mapped to the Bloch sphere. As 
is well known, any time-independent Hamiltonian acting on this system can be represented by a fixed vector, and 
the state's evolution under the Hamiltonian can simply be described by a precession of the Bloch vector for the state 
about the Hamiltonian vector. If we choose the basis to be S = {|e)|0), where the first state lies at the top of 

the Bloch sphere and the second state lies at the bottom, then we have 

^{-5/2 q \ . , r. K ^ 

where a is the (pseudo) vector associated with the Pauli matrices, and (5,0, — |) becomes the vector representing the 
Hamiltonian. Notice that all points along the "equator" of this Bloch sphere are maximally entangled, because they 
are all equivalent to ^(|e)|0) -I- by a local unitary operation. When the detuning is small, the vector for Hi 

points almost along the x-axis, and the state rotates from |e)|0) to |.g)|l) as it evolves in time, crossing the equator 
twice, so it passes through two maximally entangled states and two separable states during each rotation. Fig. [3] 
illustrates this Bloch sphere picture of the entanglement evolution, which accounts for all the qualitative features that 
appear in Fig. [51 



B. Multiple Modes 

With additional modes of the electromagnetic field, the behavior of the system quickly becomes more complex. Let 
us begin with the bipartite entanglement between the atom and the EMF when the field has only two modes. First 
consider the "symmetric" case, where gi — g2 = g and Si = —62 = S. If the initial state is pure, then we have 

5'^ + 2g^ cos (ty/2g^ + sA 

and we can easily compute the entanglement as a function of time. If the initial atomic state is mixed, we can still 
compute the LN by the technique mentioned in IIII Al The dependence on initial conditions in Fig. [4] appears quite 
similar to the single mode case, while the dependence on S shown in Fig. |4] is somewhat more complex. We are 
more interested, however, in cases where one of the field modes is very close to resonance, so a more practical case to 
investigate is that of three field modes. 

In considering a finite number of field modes, we are primarily thinking of modeling the intracavity field in cavity 
QED, particularly in the optical regime. We will adopt the simplest possible model for such a cavity, two perfectly 
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Figure 4: The dependence of the EOE on both time t and the detuning 5 from the atomic resonance frequency for the case 
where there are two EMF modes. 

conducting, infinite, parallel plates(the cavity model used in There are two classes of modes between the 

plates: the discrete set of modes with a purely longitudinal wave vector, which we will call the "cavity modes", and 
the continuum of modes with a transverse component, which propagate outward along the plates. If we ignore the 
coupHng between the atom and these transverse travehng modes but rely on the fact that the effect of a single mode 
diminishes as it becomes detuned far from resonance, then we can try to approximate the cavity by only the finite 
number of modes M . This should reflect an actual cavity in the strong-coupling limit, where g is much larger than 
the spontaneous emission rate or the cavity decay rate. 

Considering the frequencies of these cavity modes, if the mode closest to resonance is detuned by an amount 5 then 
all other modes will be detuned by an amount 5 + k/S. where A is the free spectral range of the cavity in angular 
frequency, which is A = ^ for our simple model. As noted earlier, all the gj values are proportional to each other, so 
if we call the value for the mode closest to resonance g, this sets the value for all the other coupling constants. Thus, 
we are left with only four constants whose values need be specified. 

We will consider the regime where (5 ^ A (as is often the case in cavity QED experiments); thus, one mode is 
almost in resonance, and the next nearest two modes have almost the same detuning. As a result, it will make sense 
to consider only odd total numbers of modes M in trying to understand the cavity. With this further specificity, we 
may rewrite the formula for \u[t)\. Let M — 2Q + 1, then 

\u{t)\^ 

where Wj are solutions to the equation 

(-^ - ^Su, + g-) n (^^ + A^) + 2 ^ E (l - J^f^) (- + (j^) ft + ^'^') = (^D 

We consider first the M = 3 case. When A is small, the time evolution behavior, shown in Fig. [SKa), becomes 
considerably more complex with the contributions of more frequency components to u{t). The effect of the addition 
of two more modes in the M — 5 case is shown in Fig. [HJb) . The clear intuitive interpretation that was present in 
the single mode case is not obvious in these cases. 

In optical frequency cavity QED experiments, one is generally working in the regime where ^ A g and 5 is 
smaller than or on the same order as g. In this case, other cavity modes are far enough away in frequency to have 
little effect on the atom, and the evolution for the multi-mode models becomes nearly identical to that of the single 
mode case. We can calculate the poles in this regime by writing Eq. (|2T1) in terms of the dimensionless parameters 
g/A, A/wo, and and solving it perturbatively in these small parameters. In experiments, these parameters might 
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(a) (b) 

Figure 5: The evolution with time t of the EOE between an atom in an initially pure, excited state and (a) three modes or 
(b) five modes of the EMF. To highlight the effect of additional modes, we choose the values for the detuning 5, atom-field 
coupling g, and free spectral range A so that A/g — 5, g/S = 10, and cjo/g ~ 10^, though these are unlikely to be characteristic 
of experimental values. 



typically have values on the order of 10 10 and 10 respectively. This perturbative analysis leads to the poles 

w = -^ (<5 ± V-^^ + (22) 

z. = ±.(fcA + ^) ke[l,2,...,^] (23) 

keeping up to terms linear in S/luo and terms quadratic in g/A and A/ujo (including terms of order g/ujo). This is 
the lowest order at which the poles in Eq. i(23|) yield non-vanishing contributions to \u{t)\ [23|. In the sum from Eq. 
(|20|) , the poles of the form shown in Eq. l(22l) give terms of order unity, while poles of the form shown in Eq. l(23|) 
yield terms of order g^/A^, so in this regime the effects of additional modes are small, but we have found an analytic 
expression for the behavior of the system for an arbitrary, finite number of modes. 



C. Full Intracavity Field 



Another approach toward understanding the effect of the cavity is to look at the long-time behavior of entanglement 
in the presence of the infinite field modes of the cavity. Again, we will use the simple cavity model mentioned earlier 
consisting of two infinitely large, perfectly conducting parallel plates. We can begin with a box having a square 
transverse cross-section of area A and a longitudinal length L with boundary conditions such that the field vanishes 
on the plates at the longitudinal boundaries and is periodic on the transverse boundaries. In this case we can write 



9k 



, where the dimensionless quantity A is the strength of the coupling to the overall field. As found in [jj 



A2 



27r2i 



E 



in the continuum limit where ^ ^ oo. Since this integral clearly diverges, we add an exponential cutoff by taking 
s ^ s — «e to regularize it. With this cutoff 



A2 



27r2L ^ 



n— — oc> 



'ik' s—k' e 



A2 



1 1 



-Z7r(s — ie) / L 



nLe + isl- (>,-^^{s-i()IL ■ 



Taking the Laplace transform yields approximately [T 

-zA2 \2 ,x2 



-^z In {ie'^tz) — — 



In r 



Lz\\ Lz^^fLz\ Lz 1^^^/Lz\ 



(24) 



+ Oie) (25) 



where In (r(2;)) is defined such that it has a branch cut on each ray {— n + | n e N&y G M > 0}, and 7e is the 
Euler-Mascheroni constant. Clearly this expression depends on the cutoff e, which we assume to be small. One cannot 



10 



in the end take the Hmit of e ^ for this model, even with renormahzation of the model parameters, e should be 
regarded as a phenomenological parameter that reflects the fact that at high frequencies the approximations that 
underlie our model (the rotating wave approximation, the two level approximation, etc.) must give way to other 
physics. 

Inserting this result for Jj-iz) into Eq. (fT2|) . u{t) may be obtained by finding the residues arising from the poles of 
the integrand. As found in |1J|, all poles have a negative real part. From Eq. ifTSj) one can see that the negative real 
part of each pole will cause its contribution to u{t) to die off exponentially with time, so the pole with the greatest 
real part will dominate at long times. 

If A is sufficiently small, this pole can be found by doing a perturbative expansion in A that yields 



= — zcjQ — Me(— ^wo) + O (a"') = —iil — 7 



(26) 



where luq = ^^o t- for notational convenience. Let us further define 



UJQ 



In (e'^'^ecjo) ■ 



(27) 



floe would be the value of fl in the case of the free field and, therefore, the physically observable dressed value of the 
two-level transition frequency. If we rewrite our expression for Zp in terms of this dressed atomic frequency, we may 
write 



LlIq — Hoc + 

SO then the pole from Eq. (|26l) becomes 



In {e'''enoo)+0{X'^) 



In 



+ 



■In - 



1, / Lfloc 

2^" -^ 



(28) 



+ 0(e) + 0(A4). (29) 



This perturbative solution will be valid as long as A is sufficiently small that the higher order terms can be ignored 
and no branch cut lies between —i0.oo and the pole calculated in Eq. ((29]). Clearly, the expression will not be 
valid when Lfloc/T^ is a non-negative integer, which is the condition for resonance. When floo is sufficiently close to 
resonance, we can verify numerically that there will be two closely spaced poles, as one would expect from a simple 
Jaynes-Cummings model treatment. 

In the case that the perturbative solution is valid and gives the only pole important at long times, then in that 
limit \u{t)\ (X e~^* from Eq. (|26|) . The value of 7 depends on the value of the quantities Qoo, A, e, and L. When A is 
sufficiently small, 7 may be obtained by taking the real part of Eq. l(29l) . In other cases it may be found numerically. 

Having determined the poles, if they are first order then the residue for each will simply be ^1 + ^^^^|^ ^ • Com- 
bining the expression for \u{t)\ at long times with our expressions for entanglement of a pure state yields a relatively 
simple behavior for the entanglement. The exponential decay of atomic coherence noted by [l3| leads to entanglement 
that falls off toward zero at long times. 



V. SUMMARY AND DISCUSSION 



In this work we have considered the interaction between a two-level atom and different numbers of electromagnetic 
field modes with the aim of gaining a detailed description and deeper understanding of the nature and dynamics of 
quantum entanglement between the atom and the field. We have produced exact results (given the usual atom-field 
interaction Hamiltonian derived from making the rotating wave, dipole, and two-level approximations) that go beyond 
the usual single-mode Jaynes-Cummings model and give the complete development of quantum entanglement in time 
for a system where the field is initially in the vacuum state. 

On the effect of initial conditions of the atom in the general case, we find that quantum entanglement is not 
affected by the angle of the initial atomic state on the Bloch sphere. For initially pure atomic states, we also find 
that qualitative features of the time evolution of entanglement remain the same for different values of 6. 

For the case where the cavity has only one dominant field mode, our calculations reproduce the familiar results for 
the Jaynes-Cummings model obtained before 0, H, Our result shows the expected oscillations of entanglement 
between the atom and the field, including periodic complete disentanglement in the resonant case. While many prior 
treatments assume a pure initial atom-field state or that the field mode is resonant with the 2LA, we have calculated 
the entanglement when the initial state (and, thus, the state at all times) of the atom-field system is mixed. We have 
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also shown the effects of detuning from resonance on the dynamics, giving a simple conceptual picture that accounts 
for all qualitative features. 

For the situation when the atom is coupled to a finite number of electromagnetic field modes, since our calculation 
made no weak-coupling approximation, we are able to give a detailed account of the behavior of entanglement in 
the strong-coupling regime. Indeed, if one wishes to relate this to an actual cavity QED system where losses due 
to spontaneous emission and cavity damping are present, our closed-system treatment will be useful exclusively in 
the case of strong coupling g k,7 at sufficiently early times where nt ^ 1 and 7t <C 1, wherein dissipation is 
insignificant and can be ignored. 

We have discussed how the effects of these additional modes far from resonance are small in the normal parameter 
range of current cavity QED experiments. On the other hand, if the system were placed in the super-strong-coupling 
regime, where g is comparable to the free spectral range of the cavity A, then we can see from Fig. [5] that the 
dynamics are significantly altered from those of the single-mode case. Our calculations add new, significant features 
to the behavior of the entanglement. 

It has been suggested that one may be able to reach the super-strong-coupling regime experimentally in a 
system with many ultracold atoms trapped in an optical lattice inside the cavity. Under the right circumstances, our 
calculations might apply to such a system. If the experimental realization were such that the intracavity field coupled 
almost identically to the internal states of all N atoms so that the system behaved like the Tavis-Cummings model 
fl^, then our model could still be applied if the atoms were initially prepared in the Dicke state |A^/2,— (A^— l)/2) 
[19] (i.e., a collective state with one excitation symmetrized over all atoms). In this case, the collective dipole of the 
atomic sample behaves like a single, two-level atom with enhanced coupling constant g gVN- Of course, in general 
the situation may be considerably more complex than a Tavis-Cummings type model, but at the least our results 
suggest interesting new dynamics for quantum entanglement in this super-strong-coupling regime. 

Finally, we have treated the case of a 2LA interacting with the infinite number of modes present in an intracavity 
field and given an account of the long-time behavior of their entanglement. In this model we have assumed that the 
mirrors are perfectly refiecting. There are two possible approaches to applying these results to a physical cavity with 
damping. If the damping time scale is sufficiently longer than the other important dynamical time scales then the 
damping can be ignored for times nt 1, and one can apply our results directly. Alternately, one could modify the 
treatment we have given by using cavity field modes that extend through the cavity mirrors to the exterior, as would 
exist with partially transmitting mirrors, though the entanglement calculated would still be with the entire field, not 
merely the intracavity portion. 

Though our solution for the cavity case could be evaluated numerically at any time, we pay special attention to the 
long-time behavior where an analytic result can be found. In order for this regime to be meaningfully applied to a 
physical system with k ^ 0, one requires that such long times are still significantly shorter than 1/k. The long-time 
solution we have presented requires all other contributions to u{t) to be negligibly small. These correspond to poles 
in the complex plane with negative real parts that are larger in absolute value. If the next most significant term has a 
real part —7', then the time scale that defines long times will be (7' — 7)"^. Given that we expect any other poles to 
move off to infinity in the complex plane as A ^ 0, if A is small enough the time scale that defines long times should 
become comparatively short, and there should be a window of validity where it can be considered long times and yet 
K may be ignored. 

In obtaining these exact solutions, we are able to see several features of the dynamics. With the infinite collection 
of modes, we find that the terms from Eq. lfT2|) are exponentially decaying, leading to a behavior in which the 
entanglement dissipates at long times. By contrast, when only a few field modes are included each term is purely 
oscillatory and no such decay is present. In each of the cases with an infinite collection of field modes, our exact 
calculation leads to the expression for Jl{z) that can be used to evaluate u{t) numerically in order to obtain the 
behavior of the system. To give more concrete results, we have examined the long-time behavior working out the 
result perturbatively in the coupling A. While this long-time behavior results in an exponential decay of coherence in 
the atom jl3| characteristic of Markovian behavior of the atom (now considered as an open system with the field modes 
integrated away), our exact result for the total system can provide a check on how good the various approximations 
introduced for any description of the open system could be. These features could be explored by solving the equations 
we derived here numerically in particular limits of interest. 

We addressed the strong-coupling regime above. Another important feature which our study of the full dynamics 
of the total system can provide some comparison with is the non-Markovian versus Markovian dynamics of a related 
open system. As is explained in the Introduction a la Zwanzig's projection operator paradigm, the exact and complete 
description of the total system (via an ordinary differential equation) can be transcribed to that of a distinguished 
subsystem (where one desires a detailed description) with the dynamics of all the other subsystems it interacts with 
subsumed in the non local kernels of integro-differential equation it obeys. This is the origin of non-Markovian 
dynamics with memory in a closed system. The open-system dynamics is only one step away from this, when one 
decides to coarse-grain away to varying extents the degrees of freedom of the other subsystems the distinguished 
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system interacts with. That is why the dynamics of such open systems are generically non-Markovian dynamics. The 
unitary dynamics of a closed combined system thus offers the best way to compare the relation of the non-Markovian 
dynamics of the open system with the Markovian dynamics of the same open system and assess the accuracy of 
various approximations introduced for its description. In the cases that include a finite number of field modes, one 
can obtain such subsystem dynamics quite directly from our solution for the total system dynamics. For a long time 
quantum atomic-optical experiments have operated in parameter ranges where the dynamics give relatively simple 
yet highly accurate results, so these issues have been largely ignored. But with stronger couplings, shorter times, and 
many correlated atoms present (the case of cold atoms in optical lattices is a familiar one) this luxury may not be 
available for too long. 

The ever increasing capability for exquisite control of atom-field interactions in the laboratory and the increasing 
demands of quantum information applications require ever increasing accuracy in our theoretical modeling of these 
quantum systems. In addition, optical cavity QED is now truly approaching the regime of strong coupling. By 
working out the dynamics of quantum entanglement in this simple system exactly, in the regime of strong coupling 
and without making any compromising approximations (like the popular Born-Markov approximation), we hope to 
begin to uncover some new features of this realm of physics. Ours is really just a small step in this direction. There 
remains much work to be done, both to understand the generic behaviors of these simple systems by model studies 
and to add more features to the theoretical models that can provide a closer depiction of reality, captured in the near 
future by higher precision experiments. 
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